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1 Introduction 



The diffractive deep inelastic scattering (DDIS) at HERA provide a very interesting example of 
the interplay between hard and soft aspects of QCD interactions. On one side, the virtuality of 
the photon probe is large (Q 2 3> Aq CD ), while on the other side, the scattered proton remains 
almost intact, loosing only a small fraction of its initial momentum. Its transverse momentum 
with respect to the photon-proton collision axis is also small. In addition to the scattered incident 
particles, a diffractive system forms which is well separated in rapidity from the scattered proton. 
The most important observation made at HERA is that diffractive processes in DIS are not rare, 
quite the contrary, they constitute up to 15% of deep inelastic events. What's more, the ratio of 
the diffractive and inclusive cross sections is constant as a function of energy of the ^y*p system 
or as a function of the photon virtuality. The latter fact reflects the logarithmic dependence on 
Q 2 of diffractive structure functions in the Bjorken limit. 

In the i-channel picture, the diffractive interactions can be viewed as a vacuum quantum 
number exchange between the diffractive system and the proton. In old days of Regge phe- 
nomenology such a mechanism of interactions was termed a pomeron. With the advent of 
quantum chromodynamic we gain a new way of understanding the pomeron by modelling it 
with the help of gluon exchanges projected onto the color singlet state. In the lowest approxi- 
mation, the pomeron is a two gluon exchange which is independent of energy. By considering 
radiative corrections to this process in the high energy limit, the famous BFKL pomeron [1-4] 
was constructed with a strong, power-like dependence on energy. This dependence ultimately 
violates unitarity which means that exchanges with more gluons have to be considered. A sys- 
tematic program to sum exchanges with gluon number changing vertices was formulated in [5,6] 
and developed in [7-10]. Other, somewhat more intuitive formulation, called Color Glass Con- 
densate [11-14], is based on the idea of parton saturation [15] in which deep inelastic scattering 
occurs on a dense gluonic system in the proton. In these approaches unitarization is supposed 
to change the asymptotic energy behaviour of the cross sections involving the pomeron from 
power-like to logarithmic. 

DDIS is particularly sensitive to the pomeron energy behaviour since diffractive scattering 
amplitudes are squared in diffractive cross sections. Thus, unitarization effects play more im- 
portant role than for the total cross section which is proportional to the imaginary part of the 
scattering amplitude. This observation is a basis of a successful description of the first diffractive 
data from HERA in which the diffractive system was formed by the quark-antiquark (qq) and 
quark-antiquark-gluon (qqg) systems. They can be viewed in the space of Fourier transformed 
transverse momenta as color dipoles [16, 17]. In the approach we follow in the forthcoming pre- 
sentation, the pomeron interaction is modelled by a two-gluon exchange which is subsequently 
substituted by the effective dipole-proton cross section fitted to inclusive DIS data [18,19]. In 
this way unitary is achieved. 

In an alternative approach to DDIS, the diffractive structure functions are defined in terms of 
diffractive parton distributions (DPD). They are evolved in Q 2 with the help of the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations [20-22]. Thus in the Bjorken limit, the 
diffractive structure functions depend logarithmically on Q 2 , i.e. they provide the twist-2 de- 
scription of DDIS. The theoretical justification of this approach is provided by the collinear fac- 
torisation theorem which is valid for hard diffractive scattering in ep collisions [23-27]. However, 
collinear factorisation fails in hadron-hadron scattering due to nonfactorizable soft interactions 
between incident hadrons [28,29]. Thus, unlike inclusive parton distributions, the DPD are not 
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Figure 1: Kinematic variables relevant for diffractive DIS. 



universal objects and in general can only be used for diffractive processes in the ep deep inelas- 
tic scattering. Nevertheless, the scale of nonuniversality can be estimated by applying them to 
hadronic reactions. 

The relation between the color dipole approach with the qq and qqg diffractive components 
and the DGLAP based description was studied in detail in [30]. In short, after extracting the 
twist— 2 part, the dipole approach provides (^-independent quark and gluon DPD. In addition, 
the qqg component, which was computed assuming strong ordering between transverse momenta 
of the gluon and the qq pair, gives the first step in the Q 2 -evolution of the gluon distribution. 
The twist-2 approach, which is based on the DGLAP evolution equations, extends the two 
component dipole picture by taking into account more complicated diffractive final state. In 
the performed up till now twist-2 analyses, the diffractive parton distributions are determined 
through fits to diffractive data from HERA [31] We will follow this approach with an important 
modification. 

The dipole approach teaches us one important lesson concerning the seemingly subleading 
twist-4 contribution, given by the qq pair from longitudinally polarised virtual photons (Lqq). 
Formally, it is is suppressed by a power of 1/Q 2 with respect to the leading twist-2 transverse 
contribution. However, the perturbative QCD calculation shows that for small diffractive masses, 
M 2 <C Q 2 , the longitudinal contribution dominates over the twist-2 one which tends to zero 
in this limit. The effect of the Lqq component is particularly important for the longitudinal 
diffractive structure function which is supposed to be determined from the high luminosity 
run data at HERA. Thus, we claim that it is absolutely necessary to consider the twist-4 
contribution in the determination of the diffractive parton distributions. The analysis which we 
present confirms its relevance for the prediction for , which differs significantly from that 
based on the pure DGLAP analysis. This is the main result of our paper. 

The paper is organised as follows. In Section [2] we provide basic formulae for the kinematical 
variables and quantities measured in diffractive deep inelastic scattering. We also describe the 
three contributions which we include in the description of the diffractive structure functions, i.e. 
the twist-2, twist-4 and Regge contributions. In Section [3] we describe performed fits while in 
Section 0] we present their impact on the determination of the diffractive parton distributions 
and diffractive structure functions. We finish with conclusions and outlook. 
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2 Basic formulae 



We consider diffr active deep inelastic scattering: ep — ► e'p'X, shown schematically in Fig. [TJ 
After averaging over the azimuthal angle of the scattered proton, the four-fold differential cross 
section is given in terms of the diffractive structure functions F® and F£: 

where y = Q 2 / (xbs) and s is the ep centre-of-mass energy squared. The expression in the curly 
bracket is called reduced cross section: 

Both structure functions depend on four kinematic variables (P,Q 2 ,xjp,t), defined as follows 

Q 2 + M 2 -t Q 2 

where — Q 2 is virtuality of the photon, t = (p — p') 2 < is the square of four-momentum 
transferred into the diffractive system, M is invariant diffractive mass and W is invariant energy 
of the 7*p system. The Bjorken variable xb = xjp(3. For most of the diffractive events \t\ is 
much smaller then other scales, thus it can be neglected in eqs. (jSJ)- The diffractive structure 
functions are measured in a limited range of t, thus the integrated structure functions are defined 

F£ 3 \/3,Q 2 ,x P )= dtF° L (P,Q 2 , X]P ,t), (4) 

* t'min 

The integrated reduced cross section is defined in a similar way. 

2.1 Twist— 2 contribution 

In the QCD approach based on collinear factorisation, the diffractive structure functions are 
decomposed into the leading and higher twist contributions 

= *sr } + <r 4) + • ■ • • (5) 

The twist-2 part is given in terms of the diffractive parton distributions through the standard 
collinear factorisation formulae [23,32-34]. In the next-to-leading logarithmic approximation we 
have 

F° (tw2) (x, Q 2 ,xjp, t) = S D + ^{C$ ®S D + Cf (8) G D } (6) 
F^ tw2 \x, Q 2 , x P , t) = g {Cl ®S D + CC®G D } (7) 

S G 

where a s is the strong coupling constant and C 2 V are coefficients functions known from inclusive 
DIS [35,36]. The integral convolution is performed for the longitudinal momentum fraction 

(C <g> F)((3) = f 1 dzC {15 /z) F(z) . (8) 



Notice that in the leading order, when terms proportional to a s are neglected, the longitudinal 

D( 

L 

gluon distributions, q^ D and go'- 



structure function p^ tw2 ^ = q. The functions Sp and Gd are given by diffractive quark and 



S D = Y,e}(3 {q f D (/3, Q 2 ,x 1P ,t) + q f D (p,Q 2 ,x 1 p,t)} (9) 
/=i 

G D = f3g D {l5,Q\x IP ,t) (10) 

Note that (5 = x/xjp plays the role of the Bjorken variable in DDIS. In the infinite momentum 
frame, the DPD are interpreted as conditional probabilities to find a parton with the momentum 
fraction x = (3xp in a proton under the condition that the incoming proton stays intact losing 
a small fraction xjp of its momentum. A formal definition of the diffractive parton distributions 
based on the quark and gluon twist-2 operators is given in [23,25]. 

The DPD are evolved in Q 2 by the DGLAP evolution equations [37] for which the variables 
(xjp,t) play the role of external parameters. In this analysis we assume Regge factorisation for 
these variables: 

q f D iP, Q 2 , xjp, t) = f P (x P , t) qfptf, Q 2 ) (11) 
g D (P, Q 2 ,X]p, t) = fjpixjp, t) g P {p, Q 2 ) . (12) 

For convenience, the functions q^pi/3, Q 2 ) and gjp(P, Q 2 ) are called pomeron parton distributions. 
The motivation for such a factorisation is a model of diffractive interactions with a pomeron 
exchange [38]. In this model fjp is the pomeron flux 



Fpif) l-2a P (t) 

fip(xp,t) = -j^-Xjp w , (13) 



where ajp(t) = ajp(0) + a'pt is the pomeron Regge trajectory and the formfactor 

Fl{t) = FU0)e- B ^\ (14) 

describes the pomeron coupling to the proton. We set -Ff>(0) = 54.4 GeV -2 [34], Bd = 
5.5 GeV -2 and a'p = 0.06 GeV -2 [31], while the pomeron intercept ajp(0) is fitted to data. 
The pomeron quark distributions are flavour independent, thus they are given by one function, 
a singlet quark distribution Ejp: 

qfp (p, Q 2 ) = q F iP, Q 2 ) = -L EiP (/3, Q 2 ) (15) 

where Nf is a number of active flavours. The question about Regge factorisation is an issue 
which should be tested experimentally. In our approach, the pomeron is a model of diffractive 
interactions which only provides energy dependence through the xjp-dependent pomeron flux. 
Its normalisation is only a useful convention since the normalisations of the pomeron parton 
distributions in eqs. (fTTj) are fitted to data. 

2.2 Twist-2 charm contribution 

We describe the charm quark diffractive production using twist-2 formulae for the cc pair gener- 
ation from a gluon. These are formulae analogous to the inclusive case in which the diffractive 
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Figure 2: The qq and qqg components of the diffractive system in the dipole approach. 



gluon distribution g D is substituted for the inclusive one [39]: 

F 2 L{(3,Q,x 1P ,t) = 2Pe c — / — C 2 ,l -»t^2 ]9 {z,fJ> c ,x P ,t , (16) 

27T J a/ 3 Z \Z Q z J 

where a = 1 + 4m 2 /Q 2 and the factorisation scale /x 2 = 4m 2 with the charm quark mass 
m c = 1.4 GeV. The coefficient functions read 

C 2 (z, r) = i {z 2 + (1 - z) 2 + 4z(l - 3z)r - 8z 2 r 2 } In 

+ |o {-1 + 82(1 -z) -4z(l -z)r} (17) 



C L (z, r) = -4z 2 r In — — + 2az(l - z) (18) 
1 — a 



with a = ^/l — 4rz/(l — z). The cc pair can only be produced if invariant mass of the diffractive 
system M 2 fulfils the following condition 

M 2 = Q 2 Q-l)> 4m 2 . (19) 



2.3 Twist— 4 contribution 



The computation of the twist-4 contribution, proportional to 1/Q 2 , is a nontrivial task and 
one could be tempted to assume that this contribution is suppressed at large Q 2 as in inclusive 
DIS. However, by analysing diffractive final states in the dipole approach it was found that 
for diffractive mass M 2 <C Q 2 (f3 — > 1), the twist-4 contribution dominates over the vanishing 
twist-2 one [19,40,41]. 

This observation is made on the basis of the perturbative QCD calculations in which the 
diffractive state is formed by the qq and qqg systems interacting with a proton through a colorless 
gluonic exchange which is a model of the pomeron interactions in QCD. In the simplest case, two 
gluons projected onto the color singlet state are exchanged, see Fig. [2j The computed amplitudes 
do not depend on energy in such a case which problem can be cured in a more sophisticated 
approach by modelling the dipole-proton cross section which fulfils unitarity conditions [19]. 
Independ of the details of the pomeron description, the diffractive mass (or (3) dependence 
is a genuine prediction of pQCD calculations. It appears that the leading in Q 2 behaviour 
components, qq and qqg from transverse virtual photons, vanish for (3 — > 1. This is not the case 
for the qq production from longitudinal photons (Lqq) which is formally suppressed by 1/Q 2 
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with respect to the leading components. Thus, this particular /3-dependence makes the Lqq 
contribution dominant for /3 — ► 1, see Fig. [3j 

The presence of the Lqq component has important consequence for the longitudinal diffractive 
structure function which is supposed to be determined from the HERA data. The formula given 
below is an important element in the description of F® in the region of large (3: 

o o3 r Q 2 ^-0) ,2 /n 2 

^ - isfc *- BM y>j^wL dk2 1 % P «*<*■ "> (20) 



1-0 Q 2 



where the function (fio(k,X]p) is given in terms of the dipole cross section a(xjp, r) and the Bessel 
functions Kq and Jq: 




<j} (k,xip) = k / drrK W- > J (k)a(xp,r). (21) 



Strictly speaking, eq. (|20|) contains all inverse powers of Q 2 but the part proportional to 1/Q 2 
(called twist-4) dominates. The dipole-proton cross section describes the interaction of a color 
dipole, formed by the qq or qqg systems, with a proton. Following [18] we choose 

a(xjp, r) = a {1 - exp (-r 2 Q 2 /4)} (22) 

where Q 2 = (xjp/xo) _A GeV 2 is a saturation scale which provides the energy dependence of 
the twist-4 contribution. The parameters uo = 29 mb, xq = 4 • 10 -5 and A = 0.28 are taken 
from [18] (Fit 2 with charm). This form of the dipole cross section provides successful description 
of the first HERA data on both inclusive and diffractive structure functions [18, 19]. A different 
parametrisation of a, without the saturation scale, is also given in [42-44]. We checked that a 
very similar description of FP^ was found in a recent analysis [45] based on the Color Glass 
Condensate parameterisation of the dipole scattering amplitude [46]. 

The relation between the dipole approach with three diffractive components and the DGLAP 
approach with diffractive parton distributions was analysed at length in [30] . Summarising this 
relation, the twist-2 part of the qq component gives a diffractive quark distribution. The twist-2 
part of the qqg component forms a first step of the DGLAP evolution which starts from a given 
gluon distribution. Both diffractive parton distributions do not depend on Q 2 , thus they may 
serve as initial conditions for the DGLAP equations at the scale which is not determined. From 
this perspective, the DGLAP approach offers a description of more complicated diffractive state 
with any number of partons ordered in transverse momenta. However, the pQCD calculations 
tell us that the twist-2 analysis of diffractive data should include the twist-4 contribution since 
it cannot be neglected at large (3. This is the strategy which we follow in our analysis. 

We also borrow from the dipole approach a general form in (3 of the initial quark distribution 
which vanishes at the endpoints (3 = 0, 1 (see eq. ([3T]) in which A q and C q are positive). A very 
important aspect of Regge factorisation (fTTI) can also be motivated by the dipole approach. It 
is a consequence of geometric scaling of the dipole cross section (|22p [30,47]. 

2.4 Reggeon contribution 

The diffractive data from the HI collaboration for higher values of xjp hints towards a contribu- 
tion which decreases with energy. This effect can be described by reggeon exchanges in addition 
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to the rising with energy pomeron exchange. Following [48, 49], we consider the dominant 
isoscalar (/2,w) reggeon exchanges which lead to the following contribution to F®: 

if (R) = E frt**>> *) MP, Q 2 ) ■ (23) 

R 

This contribution breaks Regge factorisation of the diffractive structure function, however, its 
presence is necessary for xp > 0.01 [50]. The reggeon flux f R is given by the formula analogous 
to eq. (USD 

f R { XF ,t) = S^e-WMI \m{t)? ^ 2aR{t) , (24) 

where aji(t) = 0.5475 + 1-4 is the reggeon trajectory. From the Regge phenomenology of 
hadronic reactions A R = 0.65 GeV and the reggeon-proton couplings are given by [49]: Fj 2 (0) = 
194 GeV -2 and F 2 (0) = 52 GeV" 2 . The functions 

\ m {t)\ 2 = 4cos 2 [7ra fl (t)/2] , \ m (t)\ 2 = 4sin 2 [vra i? X*)/2] (25) 

are signature factors for even (/jj) and odd (u) reggeons, respectively. We could also consider 
isovector reggeons (a,2,p) but their couplings to the proton are much smaller and we neglect 
them. Finally, the reggeon structure function Fr is given by [49] 

F R (J3) = A R r° m (l-P) 2 , (26) 

where the normalisation A R is a fitted parameter. Thus, in the first approximation, we neglect 
the Q 2 -dependence of the reggeon contribution. 



3 Fit details 



Collab. 


No. points 


Data 


\t -range 


Q 2 -range 


/3-range 


HI [50] 


72 


LP 


[0.08, 0.5] 


[2,50] 


[0.02,0.7] 


ZEUS [51] 


80 


LP 


[0.075 , 0.35] 


[2,100] 


[0.007,0.48] 


HI [31] 


461 


M Y < 1.6 


[l^minl j 1] 


[3.5 , 1600] 


[0.01,0.9] 


ZEUS [52] 


198 


My < 2.3 


[l^minl i Oo] 


[2.2,80] 


[0.003 , 0.975] 



Table 1: Kinematic regions of diffractive data from HERA. LP means leading proton data and 
My is invariant mass of a dissociated proton. Dimensionfull quantities are in units of 1 GeV. 

In our analysis we use diffractive data from the HI [31,50] and ZEUS [51,52] collaborations. 
In Table Q] we show their kinematic limits in which they have been measured. The minimal value 
value of \t\ is given by 

r 2 

It . I ~ F m 2 (27*) 

\ b min\ — -. d 5 K^'J 

1-Xjp 

where m p is the proton mass. The leading proton data from HI, measured in the range given 
in Table[H were corrected by the HI collaboration to the range \t m i n \ < \t\ < 1 GeV 2 . 

The ZEUS data are given for the diffractive structure function F^, thus we use in our 
analysis the following formulae 

p D = F D(t W 2) + F D(R) + f d_ (28) 

F D =F *>M + F D 9 . (29) 
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No 


Data 


Fit 


ajp(O) 


Ar 


A q 


B q 


c q 


A 9 


B 9 


c 9 


X 2 /N 


1 


HI (LP) 


tw-2 


1.098 


0.29 


1.75 


1.49 


0.5* 


2.09 


0.67 


0.80 


0.48 


2 


ZEUS (LP) 


tw-2 


1.145 


1.05 


2.13 


1.51 


0.5* 


10.0* 


1.03 


2.26 


0.40 


3 


HI 


tw-2 


1.117 


0.49 


1.33 


1.63 


0.34 


0.17 


-0.16 


-1.10 


1.04 


4 




tw-(2+4) 


1.119 


0.48 


1.62 


1.98 


0.59 


0.04 


-0.56 


-1.68 


1.17 


5 


ZEUS 


tw-2 


1.093 


0.0* 


1.68 


1.01 


0.5* 


0.49 


-0.03 


-0.40 


1.35 


6 




tw-(2+4) 


1.092 


0.0* 


1.20 


0.85 


0.57 


0.07 


-0.52 


-1.48 


1.82 



Table 2: The fit parameters to HI nd ZEUS data. The presence of twist-4 in the fits is marked 
by tw-(2+4). The parameters with an asterisk are fixed in the fits. 

The longitudinal twist-4 contribution is present on the r.h.s. of eq. (|28p since F 2 D is the sum of 
the contributions from the transverse and longitudinal polarised virtual photon. The HI data, 
however, are presented for the reduced cross section ([2]). Thus we substitute relations (|28p and 
(|29|) in there and use 

-D _ \ V D{tw2) „D(R) V 2 p D(tw2) \ 2(1 ~ V) p fl /onN 

a r~[ F 2 + F 2 - l + {1 _ y? F L | + l + {l _ y) 2 F Lm- ( 3 °) 

The expression in the curly brackets is the twist— 2 contribution while the last term is the twist-4 
one. Notice that the difference between and a® is most important for y — > 1. 

We fit the diffractive parton distributions at the initial scale Qq = 1.5 GeV 2 , assuming the 
Regge factorised form (jlip with the following pomeron parton distributions [31]: 

PS P (J3) = A q (3 B «(l-f3fi (31) 

PgrW) = A g p B »(i-pfs. (32) 

The six indicated parameters are fitted to data. We additionally multiplied both distributions 
by a factor exp{— a/(l — /?)} with a = 0.01 to secure their vanishing for (3 = 1. This factor 
is only important when C q or C g becomes negative in the fits. For the evolution, we use the 
next-to-leading order DGLAP equations with Aqcd = 407 MeV for Nf = 3 flavours [53]. 

The pomeron flux in eq. (Ilip is integrated over t in the limits given in Table Q] which leads 
to the form 

MX1P) = SI { e ~ Bltm " A ~ e " B|W1 } 4" 2QJP(0) • (33) 
The shrinkage parameter B equals 

B = B D + 2a' ]P \n{l/x I p) (34) 

with B D = 5.5GeV~ 2 and a' F = 0.06 GeV~ 2 [50]. 

In summary, we have eight fit parameters altogether: the pomeron intercept ajp(0), reggeon 
normalisation Ar in eq. (|2"B"P and six parameters in eqs. (|31l32p 

4 Fit results 

The data sets from Table[T]were obtained in different kinematical regions, using different methods 
of their analysis. Thus, we decided to perform fits to each data set separately. The values of 
the fit parameter are shown in Table [2j The difference between them can be attributed to the 
scale of uncertainty of our analysis. In each case we preformed two fits: with and without the 
twist-4 formula added to the twist-2 contribution. 
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4.1 Leading proton data 



We started from fits to the leading proton data. The fit parameters in this case are displayed 
in the first two rows of Table El We only show the twist-2 fit results since they are not changed 
in fits with the twist— 4 term. This happenes because the leading proton data comes from the 
region of (3 values where the twist-4 contribution is small (/? < 0.7 for HI and (3 < 0.5 for 
ZEUS), see Fig. [3 

The data with a dissociated proton (DP) which are measured in the region of large (3 influence 
most the value of the parameter C g which controls the behaviour of the gluon distribution at 
P — ► 1. For the LP data C g is positive and the gluon distribution is suppressed near (3 ph 1, while 
for the DP data C g is negative and the gluon distribution is strongly enhanced. This shows that 
the data with (3 > 0.7 are crucial for the proper analysis. Without this kinematic region we lose 
important information about diffractive interactions. Thus, from now on we concentrate on the 
analysis of the DP data. 

4.2 HI data 

The fit parameters to the HI data with a dissociate proton are given in the third and fourth rows 
of Table [21 We see that the fit quality is practically the same for both fits, with and without the 
twist-4 contribution. The presence of the reggeon term improves fit quality by 30 units of x 2 
for 461 experimental points. A good quality of the fits is illustrated in Fig. [J] which also shows 
that the reduced cross sections (|30j) from the twist-2 (solid lines) and twist-(2+4) fits (dashed 
lines) are very close to each other. 

In Fig. [5] we show our results on the reduced cross section for the largest measured value of 
(3 = 0.9. In this region, the twist-4 contribution, shown by the dotted lines, cannot be neglected. 
We see that the curves from both the twist-2 (solid) and twist-(2+4) (dashed) fits describe data 
reasonable well. However, the curves with twist-4 have a steeper dependence on xjp (energy) 
than in the pure twist-2 analysis. This observation is by far more pronounced in the analysis of 
the ZEUS data performed for the structure function . 

The diffractive parton distributions from our fits are shown in Fig. [6] in terms of the pomeron 
parton distributions, /?Ejp(/3, Q 2 ) and (3gjp((3, Q 2 ). Being independent of the pomeron flux, such 
a presentation allows for a direct comparison of the results from fits to different data sets. We see 
that the singlet quark distributions are quite similar while the gluon distributions are different. 
In the fit with twist-4, the gluon distribution is stronger peaked near (3 ~ 1. This somewhat 
surprising result can be understood by looking at the logarithmic slope of F® for fixed values 
of (3. From the LO DGLAP equations we have schematically: 

at? ~ ataff = + p » c8G - - W p « < 35 > 

where the negative term describes virtual corrections. For large (3, the measured slope is negative 
which means that the virtual emission term must dominate over the real emission ones. The 
addition of the twist-4 contribution to F^, proportional to 1/Q 2 , contributes a negative value 
to the slope which has to be compensated by a larger gluon distribution in order to describe the 
same data. 

In Fig.[7]we present our most important results. On the left panel, the F^ 1 structure function 
is shown from both fits, with and without the twist-4 contribution (shown by the dotted lines). 
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We see no significant difference between these two results. However, the longitudinal structure 
function differs significantly for the two fits (right panel) due to the twist-4 contribution. 
Let us emphasise that both sets of curves were found in the fits which well describe the existing 
data on a® , including the large (3 region. Thus, an independent measurement of in this 
region would be an important test of the QCD mechanism of diffraction. 

4.3 ZEUS data 

The results of same fits performed for the ZEUS data are shown in the last two rows of Table [2] 
This time the Regge contribution (|26|) is not necessary since fits give the reggeon normalisation 
Ar ~ 0. In general, the fit quality is worse than for the HI data. 

As shown in Fig. El the biggest difference between the twist-2 and twist-(2+4) results occurs 
at large (3 values. This is analysed in detail in Fig. [9] We see that the presence of the twist-4 
term in the fit (dashed lines) improves the agreement with the data in this region. In particular, 
a steep dependence of F® on xjp is better reproduced by the twist-(2+4) fit then by the twist-2 
one (solid lines). This dependence is to large extend driven by the twist-4 contribution (dotted 
lines). 

The behaviour of the diffractive parton distributions and structure functions, shown in 
Figs. [10] and [XT] respectively, is very similar to that found for the HI data. The gluon dis- 
tribution from the fit with twist-4 is stronger peaked near (3 ~ 1 and the longitudinal structure 
functions in the large f3 region is dominated by the twist-4 contribution. 

We summarise the effect of the twist-4 contribution in Fig. [12] showing the predictions for 
the longitudinal diffractive structure function F^ . Ignoring this contribution, we find the two 
solid curves coming from the pure twist-2 analysis of the HI (upper) and ZEUS (lower) data. 
With twist-4, the dashed curves are found, the upper one from the HI data and the lower one 
from the ZEUS data. There is a significant difference between these two predictions in the 
region of large j3. We believe that the effect of the twist-4 contribution will be confirmed by the 
forthcoming analysis of the HERA data. 

5 Conclusions 

We performed fits of the diffractive parton distributions to new diffractive data from the HI and 
ZEUS collaborations at HERA. In addition to the standard twist-2 formulae, we also considered 
the twist-4 contribution which dominates in the region of large (3. This contribution comes from 
the diffractive production of the qq pair by the longitudinally polarised virtual photons. The 
effect of the twist-4 contribution on the diffractive parton distributions and structure functions 
was carefully examined. The twist-4 contribution leads to the gluon distribution which is peaked 
stronger at [3 ~ 1 than in the case without twist-4. 

The main result of our analysis is a new prediction for the longitudinal diffractive struc- 
ture function F£. The twist-4 term significantly enhances F^ in the region of large [3. A 
measurement of this function at HERA in the region of large (3 should confirm the presented 
expectations which are based on the perturbative QCD calculations. The obtained diffractive 
parton distributions can also be used in the analysis of diffractive processes at the LHC, in 
particular, to the estimation of the background to the diffractive Higgs production, see [54] for 
a recent discussion. 
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Figure 3: Three contributions to from: qq and qqg from transverse (T) and longitudinal (L) 
photons [19] for xjp = 0.003. The twist-4 contribution Lqq is indicated by the yellow band. Old 
ZEUS data points are shown. 
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Figure 4: Reduced cross section a r for HI data as a function of xjp. Solid lines: twist— 2 fit, 
dashed lines: twist-(2+4) fit. 
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Figure 8: Diffractive structure function F 2 as a function xjp for ZEUS data. Solid lines: 
twist— 2 fit, dashed lines: twist— (2+4) fit. 
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Figure 12: Predictions for (3) for x F = 1CT 3 and Q 2 = 10 GeV 2 from the twist-(2+4) fits to 
the HI (upper dashed line) and ZEUS (lower dashed line) data. The solid lines show predictions 
from pure twist-2 fits to the HI (upper) and ZEUS (lower) data. 
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